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The correlation function factorization 



(4>(xx)4>(x2)lp{z, Z)) = K^/{4>{x ± )^(z,z)){(t){x2)'4'{z,z)){(t){x 1 )<j){x2)) , 

with K a boundary operator product expansion coefficient, is known to hold for certain scaling 
operators at the two-dimensional percolation point and in a few other cases. Here the correlation 
functions are evaluated in the upper half-plane (or any conformally equivalent region) with xi and 
X2 arbitrary points on the real axis, and z an arbitrary point in the interior. 

This type of result is of interest because it is both exact and universal, relates higher-order 
correlation functions to lower-order ones, and has a simple interpretation in terms of cluster or loop 
probabilities in several statistical models. 

This motivated us to use the techniques of conformal field theory to determine the general condi- 
tions for its validity. 

Here, we discover that either (</>i, 3(xi)<f>i, 3{x2)4>i/2. o(z, z)) or {4>3,i(xi)4>3,i(x2)4>o,i/2(z, z)) factor- 
izes in this way for any central charge c, generalizing previous results. In particular, the factorization 
holds for either FK (Fortuin-Kasteleyn) or spin clusters in the Q-state Potts models; it also applies 
to either the dense or dilute phases of the O(n) loop models. 

Further, only one other non-trivial set of highest-weight operators (in an irreducible Verma mod- 
ule) factorizes in this way. In this case the operators have negative dimension (for c < 1) and do 
not seem to have a physical realization. 

Keywords: correlation functions, factorization, percolation 



I. INTRODUCTION 

Usefully characterizing fluids at thermal equilibrium requires knowledge of correlation functions. Many expressions 
for important experimental and theoretical quantities involve these functions. Higher-order correlations of quantities 
such as the density p(x) at several points, e.g. {p(x\)p(x2)p{xz)...) (where the brackets denote a thermal average) are 
ubiquitous but especially difficult to calculate. A variety of approaches have been proposed (for a review, see [1]). 
One idea is to factorize the higher-order correlations in terms of lower-order correlations (with fewer points), since 
the latter are generally better known. 

Recent work, using conformal field theory, exhibited several exact formulas in which three-point correlations factorize 
in terms of two-point correlations or correlations involving one point and an interval [2, 3]. Related results (involving 
one point and two intervals) give factorizations that are not quite exact, but with small corrections [4], To our 
knowledge, these are the only such exact (or almost exact) results for interacting fluids of any type. These formulas 
apply to percolation in two dimensions at the percolation point in the upper half-plane (or any simply connected 
5-j \ region) [2, 3, 4], the critical Q-state Potts models [5], and critical percolation in three dimensions [6]. 

One of these exact factorizations, first found in [2] for a system exhibiting critical percolation, may be understood in 
terms of the probability that various points are connected via clusters, as shown in Figure 1 (for a cluster connecting 
all three points). Here 



V(xi,X2,z) = Ky/V(xi,X2)P(xi,z)'P(x2,z) , (1) 
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FIG. 1: (color online) Cluster connections for factorization. 

where x\ and X2 are distinct boundary points (e.g. on the real axis), z lies on the boundary or in the interior, and V(-) 
is the probability that the given points lie in the same cluster. The result holds in any region conformally equivalent 
to the upper half-plane, with K a universal constant. 

The factorization in (1) superficially resembles the Kirkwood superposition approximation [8]. However, the for- 
mulas are not the same, and (1) is both exact and universal, while there does not seem to be any interacting fluid for 
which the Kirkwood formula applies exactly. 

The result (1), for critical percolation, corresponds to a c = CFT correlation function with boundary and bulk 
magnetization operators = 13 [9] and tp = 0i/2,o [10]. In this case the constant K has been evaluated [3, 7], giving 

2 7/2 n 5/2 

K ~ 3 3 / 4 r(i/3) 9 / 2 

= 1.0299268... , (2) 

and both (1) and (2) agree with high-precision numerical simulations [2]. Other factorizations are also known to agree 
well with numerical results [2, 4, 5, 6]. 

In terms of conformal quantities, (1) can be expressed as 

(cP( Xl )4,(x 2 )ij(z, z)) = Ky/WfaWz,*)) (</>(x 2 )ip(z,z)) Wrttf&f) , (3) 

where and tp are appropriately chosen boundary and bulk primary operators, respectively, and K is a (boundary) 
operator product expansion coefficient. 

In this paper, we use conformal field theory to find all choices of and ip for which the factorization (3) applies. We 
show that there are only three non-trivial combinations of highest-weight operators (of irreducible Verma modules) 
that factorize in this way. In each case (3) holds for general central charge. Operators that can be identified for 
continuous central charge are given Kac indices that do not change with c, as is common for continuous parameter 
loop models. Thus the bulk magnetization operator mentioned above is denoted 0i/2.oi even though neither index is 
a positive integer. 

The percolation model is equivalent to the FK random cluster representation of the Q = 1 Potts model. The 
equivalent correlation functions for the < Q < 4 Potts models factorize in the same way [5] . Explicit values of K 
are known for these cases [7], though the probabilistic interpretation similar to (1) is slightly more complicated. Thus 
this factorization extends to FK clusters in Potts models with < Q < 4. 

The factorization also applies when = 03 : i and ip = 0o.i/2- These weights are Potts model exponents at the 
tricritical point. There is a tricritical Q-state Potts model that encodes information about the geometrical clusters 
of the Potts model, but at a different Q value [11] (tricritical Q = 1 corresponds to geometric Q = 2, for example). 
Thus this second solution implies that the geometric clusters (blocks of like-spin in the Q £ Z case) obey the same 
factorization as the FK clusters. From the Schramm-Loewner Evolution (SLE) point of view, the FK cluster results 
are dual to the geometric cluster results under k — > 16/ n, where k is the SLE parameter. Using the above, we see 
that the previous result [2] (at k = 6) extends to 8/3 < k < 8. Note that the connection between FK and geometric 
clusters is identical to that between the dense and dilute phases of the 0(n) model. 

Finally, the factorization holds for cf> = 01/3.1/3 and ip — 02/3,2/3- Here the dimensions of both operators are 
negative for c < 1. These operators have a simple explanation in the vertex operator formalism, but we are not aware 
of any physical realization. 
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II. CONDITIONS ON THE CONFORMAL WEIGHTS 

In this section we make use of the operator product expansion to determine necessary conditions on the conformal 
weights for factorization. We then solve the conditions, which is sufficient to determine all cases for which (3) is valid. 
The operators in (3) satisfy 

(j){x)(j){x + Ax) = C Hq Ax h «~ 2h * £ Ax^/3 {e} L_ {e} q(x) , and (4) 
ges {£} 

^{z, z) = y h "~ 2hi - E y W P{i}L- { t } p(x) , (5) 

PER {£} 

where S is the set of all operators appearing in the boundary-boundary fusion rules, and R is the set of all operators 
appearing in the bulk-boundary fusion rules. We assume that the boundary fields q and p on the right hand side 
of these relations are elements of standard irreducible Verma modules so that the normal scheme of generating 
descendants via the Virasoro generators applies. Note that this assumption excludes logarithmic modules. However, 

1/2 

if the operators were logarithmic, the right hand side of (3) would generate log ' (•) factors, which are incompatible 
with the general form of a logarithmic four-point function. 

Now the orthogonality of conformal families ensures that the only bulk-boundary fusion channel that leads to 
non-zero three-point functions (<j)(xi)i/>(z, z)) on the right-hand side of (3) is the <fi channel. Similarly, only channels 
common to both (4) and (5) will contribute to the left hand side of (3); i.e. only channels in Sfl R will contribute to 
the four-point function (^)(xi)cf)(x2)ip{z, z)). 

Thus a necessary condition for the factorization is <p € S fl R. Hence only the </> channel remains on the right hand 
side, and if one expands it about y = 0, this determines the powers of y that enter. It follows that only that channel 
can contribute to the left hand side (otherwise different powers of y would appear). Thus, we have <f> = S D R, i.e. 
only the 4> channel can be common to the two operator product expansions. This condition can also be expressed in 
terms of operator product expansion coefficients as 

<f> = Q n P «-» C^pC^. p tx 5(p, 0) . (6) 

This is a strong condition on the permissible fusion rules, as we will see. 
Limiting the calculation to the <f> channel only gives 

4>{x)4>{x + Ax) = C u<t> Ax'' 1 * Ax W f3{i}L- {e} <j){x) , (7) 

TO 

which we insert into the left hand side of the factorization (3) giving 

(<f>( Xl )iKx2)il>(z,z)) = C^ U Ax- h *Y^ xW Pw CX -W^ x ^^~ z )) . ( 8 ) 

TO 



where 



v / (i-i)fr d Zi \ = (e-i)hi (£-i)hj a s 

- l ^\{zi-xY {z % -xY' 1 J {z-x) 1 (z-x) 1 - 1 (z-x) 1 (z-xY- 1 (> 

is the standard operator arising from commuting L_ /« over all the fields in the correlator and Ax := x<i — X\. 
The form of the two point function in expression (8) is fixed by conformal symmetry 



(#*i)V>M> = CWM — ) { {z _ Xl)i z_ Xl) ) , (io) 



where encodes the boundary operator normalization [12]. Inserting the explicit form of (9) in (8) and making use 
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of (10) gives, for the left hand side of the (3), 



^(x{]^(x 2 )^(z,z)) = C 000 Ax- 11 * (ffaMz, z)) l + Ai/3; 



2 d^jx^jz^)) C x J 2 (<j)(xi)^(z,z)) 
Cm Ax~ h * (4>(xi)1>{z, z)) (1 + Ax pxh^ (—^— + ^-—) ■ ■ ■ (11) 

\ \Z — X\ Z — X\J 



V(z-xi) 2 (z-xi) 2 (z - a;i)(z - Xl) 

+ Ax 2 p 2 ( + + £*±±± + h *~^ ) + (Axf 
\(z-xiy [z-xxY (z — X\)(Z — X\) ) 

where C x _ 1 = d x was used to simplify the expressions. 

Now consider the right hand side of (3). One of the correlation functions appearing there is given in (10) and the 
others are 



(<j>(xi)^(x 2 )) = ^Ax~ 2h * , (12) 
{<f>(x2)ip(z, z)) = C^a^ 



z — z\ ,l * -2 ' 1 "'' / 



2i J \(z - x 2 )(z - x 2 ) 

z — z\ v w 1 



2i J \(z — xi — Ax)(z — xi — Ax) 

Ax \ ~ h * ( Ax ' 



= (^^^ll- — I (13) 
\ z-xij \ z-xxj 

Combining all three gives an exact expression for the right hand side of (3): 



\ z—Xi J \ z—xi J 



Comparing the terms of O(Ax ) in (11) and (14) we find that K = C^^a^ 1 / 2 . Comparing at ©(Ax 1 ) and 0{Ax 2 ) 
lead, respectively, to 

(-*— + = pih^ (-L- + ^— ) (15) 

2 \Z — Xl Z — Xl J \z — Xl Z — Xl J 



\ (z — xi) 2 (z — xi) 2 (z — Xi)(z — xi) J \ (z — xi) 2 (z — xi) 2 (z — xi)(z — xi) 

h,/, + hjf, + hfj, — 2h,p 



+/?: 



(z — xi) 2 (z — xi) 2 (z — xi)(z — xi) 



giving 



1/2 = Pi , (17) 

114(2 + hf) = 8/3 n h^(l + hj,) + 8/32 (/ty + M > and ( 18 ) 

V = 8/3n^ + 4/3 2 (/i -2^) , (19) 

as necessary conditions for the factorization. 

On the other hand, we can directly determine the beta coefficients in (7) using the methods of [13], giving 

1/2 = Pi , (20) 

(1 + V) = 4(1 + 2fy)/9 u + 6j3 2 and, (21) 

Ah^ = Uh^u + (8h<t, + c)/3 2 , (22) 
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where c is the central charge. This analysis is valid as long as </> is a regular operator whose descendants have a 
standard irreducible Verma module structure. This assumption was already made above, so there is no additional 
loss of generality. 

(21) and (22) may be inverted as long as h$ ^ hip, ^2,1 • (The condition reflects the fact that these two operators 
have second order null vectors.) However for these operators [<fi] ^ [(/>] x [<fi] so they are not of interest to us. For any 
other case we find 

a c(l + V) + 8/i (^ - 2) 

Pn 4c(l + 2h+) +8^(8^ -5) ' 1 j 

« = M 1 + 5 M f24) 

c(l + 2h+) + 2fy(8fy - 5) ' 1 ' 

Inserting these expressions for /3n and /?2 into (18) gives either 

^ = , (25) 

or 

- 8(1 + 5M ■ (26) 

(Note that (25) or (26) satisfies (19) as well.) Either one is a necessary (but not sufficient) condition for the factorization 
(3) to hold. Note that if, following (25), we select the identity for cj> the factorization becomes trivial. 

Continuing the above analysis through terms of 0(Ax 3 ), we find no new conditions. At 0(Ax 4 ), three equations 
emerge. Since they are rather long we refrain from displaying them. However, making use of (26) in each of them in 
order to re-express in terms of h^ gives rise to the single new condition 



hi(3ht + 2)(c - 27/i - 1) ((c - 7)fy + c + 3/4 + 2) 
(5/i + l)(c + 8/i - l)(5c(2/ty + 3) + 2(/i - l)(8/i - 33)) 



. (27) 



There are five solutions of (27), but the apparent solution = —2/3 fails if we continue the analysis to 0(Ax 6 ). 
We next examine the remaining four solutions. 

First, hs = 0, which allows any value for h^. This is the trivial solution with <p equal to the identity. It obeys the 
factorization with K = 1. 

Next is 2 + c — (7 — c)h + 3h 2 = 0, solved by the dimensions of the Kac operators ^1.3 and 03,1. By (26) 
the corresponding bulk operators have dimensions = ft-o,i/2 and ft.1/2.0 respectively. As we mentioned in the 
introduction these two pairs h$ and are critical and tricritical Potts model exponents. 

One might question these solutions, since the analysis at 0(Ax 4 ) is singular for both <fi = ^1.3 and 03. 1, due to 
level three null descendants in these conformal families. However, it is straightforward to explicitly construct the 
differential equation associated to the third order null state and verify the factorization directly (as in [2] and [5]). 

These two correlation functions generalize the results of [2], which considered percolation only. The appearance 
of both critical and tricritical exponents means that the factorization (3) holds for either FK clusters or geometric 
clusters in the Q-state Potts models for < Q < 4. Note that in the former case, the boundary conditions are free 
(except at the points X\ and X2), while for the latter they are more complicated (the 3-state Potts model, for example, 
is discussed in [14]). 

The factorization also has an interpretation in the 0(n) loop model. Here, the geometric magnetic correlation 
between two boundary and one bulk points is measured, i.e. whether there is a path connecting all three points that 
does not cross an 0(n) loop. The two sets of correlation functions for each central charge represent the dense phase 
(corresponding to the FK clusters) or the dilute phase (corresponding to the geometric spin clusters). For a given 
central charge these two phases are dual in the SLE sense. The factorization holding for dual sets of operators (where 
the order of all Kac indices in the correlation function are reversed, or equivalently, the SLE parameter k — > 16/re) 
therefore implies that it applies across the SLE or conformal loop ensemble transition from dilute to dense loops that 
occurs at k = 4. In this sense these two sets of Kac operators constitute a single physical factorization valid for all 
8/3 < k < 8. 

The remaining solution contains the boundary dimension, = —(1 — c)/27 = /ia/3,1/3 which implies the bulk 
conformal weight, = —5(1 — c)/216 = /i2/3,2/3- Both of these are negative for c < 1. 

These weights are unique as Coulomb gas vertex operators because they are the only weights which satisfy the 
conditions for this factorization at arbitrary central charge without screening operators. By contrast, the case discussed 
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above with h$ = /11.3 and hj> = hi/2,0 requires two screening operators. The absence of screening charges means that 
this collection of operators will satisfy the factorization in question exactly for arbitrary central charge. However, it is 
not immediately apparent what these operators may represent in terms of cluster or loop degrees of freedom through 
Coulomb gas or SLE formalisms. 

Finally, note that the analysis performed to 0(Ax 4 ) restricts the boundary conformal weight to one of five values, 
of which four (including the trivial case of the identity operator) survive to higher order. This manipulation is valid 
only for conformal families with standard Verma module structure up to and including level four. Therefore Kac 
operators with null states at levels one through four could conceivably factorize without satisfying the condition (27). 
However, there are only eight operators of this type, and of these only fii.i, 0i j3 , and 03,1, already considered, obey 
the fusion constraint [(f)] x [(f)] = [(f)] . Hence no new sets of operators arise. 

III. CONCLUSIONS AND DISCUSSION 

In this note, in the context of conformal field theory, we investigate the uniqueness of the factorization (3) origi- 
nally discovered in the case of critical percolation. We restrict our attention to highest weight operators in a standard 
Verma module. Thus logarithmic modules are not explicitly considered (however, we do not expect correlation func- 
tions involving logarithmic operators to be consistent with (3) in general). Using techniques of conformal field theory, 
we show that only three non-trivial correlation functions of conformal operators can factorize in this way. They are 

(<t>iA x i) < f>iA x 2)<i>i/2,o( z 'Z))> (h,i( x i) ( l ) 3,i{x2)(f>o,i/2(z,z)), and (</>i/ 3 ,i/3(a ; i) ( / ) i/3,i/3( a; 2)02/3,2/3( z ^)>- These factor- 
izations hold for arbitrary central charge. The first two apply to FK or spin cluster connection probabilities in the 
Q-state Potts models, respectively, and to both phases of the 0(n) model. The last one has negative weights for c < 1 
and, to our knowledge, no physical realization. 

Note that one important condition for factorization is the existence of algebraic solutions for the correlation func- 
tions. Such solutions occur in conformal field theory in various cases. Some of these give rise to factorizations different 
from (3), see for instance [3]. It is probable that some of these can be generalized as well. However, in this work we 
have restricted our attention to finding all instances of (3), since to our knowledge it is the simplest exact factorization. 

Of course, our results do not address the interesting question of why the factorization is valid, or whether there is 
a single explanation that holds for the various different physical models to which it applies. 

It may prove interesting to consider factorization in higher dimensions, but the machinery of CFT is considerably 
weaker in dimensions greater than two and different approaches would be necessary. We plan to explore factorization of 
the bulk-boundary-boundary three point correlation in critical percolation in three dimensions, along with numerical 
simulations, in future work [6]. 
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